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ABSTRACT 

For the first time it has become possible to compare global 3D general relativistic magnetohydrodynamic 
(GRMHD) jet formation simulations directly to very-long baseline interferometric multi-frequency polarization 
observations of the pc-scale structure of active galactic nucleus (AGN) jets. Unlike the jet emission, which re- 
quires post hoc modeling of the non-thermal electrons, the Faraday rotation measures (RMs) depend primarily 
upon simulated quantities and thus provide a robust way in which to confront simulations with observations. We 
compute RM distributions of 3D global GRMHD jet formation simulations, with which we explore the depen- 
dence upon model and observational parameters, emphasizing the signatures of structures generic to the theory of 
MHD jets. With typical parameters, we find that it is possible to reproduce the observed magnitudes and many 
of the structures found in AGN jet RMs, including the presence of transverse RM gradients. In our simulations 
the RMs are generated within a smooth extension of the jet itself, containing ordered toroidally dominated mag- 
netic fields. This results in a particular bilateral morphology that is unlikely to arise due to Faraday rotation in 
distant foreground clouds. However, critical to efforts to probe the Faraday screen will be resolving the transverse 
jet structure. Therefore, the RMs of radio cores may not be reliable indicators of the properties of the rotating 
medium. Finally, we are able to constrain the particle content of the jet, finding that at pc-scales AGN jets are 
electromagnetically dominated, with roughly 2% of the comoving energy in nonthermal leptons and much less in 
baryons. 

Subject headings: galaxies: jets - magnetic fields - polarization - radiative transfer - radio continuum: general - 
techniques: polarimetric 
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1. INTRODUCTION 

Spatially resolved radio images of the polarized emission 
of active galactic nucleus (AGN) jets on pc-scales, obtained 
via very-long baseline interferometry (VLBI), have provided 
a wealth of information regarding the large-scale structure of 
jet magnetic fields within the emission region. For many ob- 
jects these observations suggest the presence of large-scale, 
radially-correlated, toroidally dominant magnetic fields within 
the emission region (Udomprasert et al. 1997; Lister et al. 
1998; Gabuzda et al. 2000; Pushkarev et al. 2005; Jorstad et al. 
2005), in broad qualitative agreement with simple analytical 
models of relativistic magnetized outflows (e.g., Laing & Bri- 
dle 2002b,a; Lyutikov et al. 2005; Zakamska et al. 2008). 

Performing these observations often requires removing the 
rotation of the polarization angle, 'I', associated with Faraday 
rotation, a consequence of propagation through intervening, 
non-radiating plasma. If the entirety of the Faraday rotation 
occurs in the foreground, ignoring finite beam effects, the ob- 
served polarization angle at a wavelength. A, is related to the 
intrinsic polarization angle via. 



fnh. = *int + RMA-. 



(1) 



where the Faraday rotation measure, RM, is independent of 
wavelength and has a value that depends upon the plasma pa- 
rameters along a given line of sight. In practice, eliminating 
the Faraday rotation requires imaging the source at a number of 
wavelengths, accurately registering each image to an absolute 
reference frame, verifying at each independent position within 
the image that the polarization angles follow the expected 
law and fitting to obtain a map of the RMs across the image. 



Despite the difficulty in performing these measurements, pc- 
scale RM maps for a number of blazars and BL Lac objects can 
now be found in the literature (Reynolds et al. 2001; Zhang & 
Nan 2001; Zavala & Taylor 2001; Asada et al. 2002; Zavala & 
Taylor 2002, 2003; Gabuzda & Chernetskii 2003; Cotton et al. 
2003; Zavala & Taylor 2004; Asada et al. 2004; Gabuzda et al. 
2004; Zavala & Taylor 2005; Asada et al. 2008a,b,c; Gomez 
et al. 2008; O'Sullivan & Gabuzda 2009b; Kharb et al. 2009; 
Croke et al. 2010). 

Typical values of the RMs of AGN jets range from roughly 
Iradm"^ to 10'*radm"^, with extreme core RMs of lO^radm"^ 
being observed in a handful of sources (though these are gener- 
ally at the centers of cooling core clusters). While the measur- 
ing of RMs allows the elimination of a foreground effect upon 
the structure of the polarization of AGN jets, it also presents 
an opportunity to study the intervening, non-radiating plasma. 
The source of the RMs is not a priori clear, with potential con- 
tributions arising from propagation through the Galactic inter- 
stellar medium, intracluster medium, the galactic halo of the 
AGN, the AGN broad-line region and through the circum-jet 
environment. It is this last possibility, propagation through the 
plasma immediately surrounding the jet, that is the subject of 
this paper. 

Free-free absorption from pc-scale foreground clouds has 
been detected in a handful of sources, implying that foreground 
clouds exist and may be responsible for a significant fraction 
of the observed RMs in some systems (e.g., Cen A and NGC, 
Jones et al. 1996; Walker et al. 2000), though there is no rea- 
son to expect these to result in ordered RM maps exhibiting 
gradients that are correlated with other aspects of the jet struc- 
ture. In contrast, there are a number of reasons to think that the 
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circum-jet environment is responsible for substantial fraction 
of the AGN jet RMs. Firstly, there is considerable structure 
on the scale of the resolved jet, which is unlikely to arise from 
Galactic or intergalactic propagation. Second, in many sources 
there is a clear gradient in the RM values along the jet axis, 
steeply decreasing with increasing distance from the radio core. 
Finally, in a handful of objects (e.g., 3C 273, Zavala & Taylor 
2001; Asada et al. 2002; Gabuzda et al. 2004; Zavala & Taylor 
2005; Asada et al. 2008a,b) there is now strong evidence for 
transverse RM gradients, which have thus far been interpreted 
in the context of a circum-jet Faraday screen containing large- 
scale toroidal magnetic fields (Blandford 1993; Gabuzda et al. 
2000; Broderick & Loeb 2009). 

Magnetic tower models are often invoked to provide a rea- 
sonable basis for the presence of helical fields in relativistic 
jets. However, models identified as magnetic tower types in 
the literature are non-relativistic models (Koenigl & Choud- 
huri 1985; Lynden-Bell 1996; Nakamura et al. 2001; Lovelace 
et al. 2002) and thus are not directly applicable to astrophys- 
ical jets (Lyutikov et al. 2005). Nevertheless, non-relativistic 
simulations of magnetic tower models suggest that the RMs 
are strongly influenced by the toroidal magnetic field, and 
correlated with the large-scale jet structure (Nakamura et al. 
2001; Kigure et al. 2004; Uchida et al. 2004; Lapenta & Kro- 
nberg 2005a; Nakamura et al. 2008). Comparisons between 
idealized, qualitative magnetohydrodynamic (MHD) models 
(Laing et al. 2006b, a) and observations of kpc-scale polarized 
emission from FRI jets suggest that in addition to an ordered 
toroidal magnetic field, at these distances such jets must also 
contain a smaller, though still significant, disordered poloidal 
field. However, at such large distances FRI jets are strongly af- 
fected by their interaction with the external medium, and thus 
these results may have little to do with the pc-scale, intrinsic 
structure of AGN jets in general. 

Self-consistent quantitative models of idealized relativistic 
MHD jets have been obtained more recently using a com- 
bination of analytical modelling (Vlahakis & Konigl 2004) 
and time-dependent simulations (Komissarov et al. 2007; 
Tchekhovskoy et al. 2009a,b). However, such idealized sim- 
ulations cannot consider the self-consistent formation of the 
jet by the black hole interacting with an accretion disk or the 
effects of magnetic field geometry and turbulence within the 
disk. Idealized studies cannot, e.g., self-consistently determine 
the effects of black hole spin on the jet power for a given mass 
accretion rate. In contrast, modern general relativistic mag- 
netohydrodynamical (GRMHD) simulations of accreting black 
holes allow a study of the self-consistent generation of rel- 
ativistic jets from rotating, accreting black holes (McKinney 
2006; Hawley & Krolik 2006; Beckwith et al. 2008; McKin- 
ney & Blandford 2009). GRMHD simulations allow one to 
study the conditions under which relativistic jets are produced, 
including the role of the magnetic field geometry near the black 
hole and how MHD turbulence affects the formation of the jet 
from the black hole or the wind from the turbulent disk (McK- 
inney & Gammie 2004; Beckwith et al. 2008; McKinney & 
Blandford 2009). Such simulations find that an ordered (but not 
necessarily large-scale) dipole-like magnetic field within the 
accretion disk leads to a quasi-steady relativistic outflow pow- 
ered by the Blandford-Znajek mechanism (Blandford & Znajek 
1977; McKinney & Gammie 2004). While such jets are dom- 
inated by a toroidal magnetic field, they remain stable against 
the disruptive helical kink (screw) mode due to relativistic rota- 
tion of magnetic field lines, sideways (lateral) expansion, finite 
mass-loading, and non-linear saturation (McKinney & Bland- 



ford 2009). 

Synchrotron emission maps have been produced for highly- 
idealized relativistic magnetized jet models (Lyutikov et al. 
2005; Zakamska et al. 2008; Gracia et al. 2009), but they gener- 
ally have major simplifying assumptions, such as perfect cylin- 
drical geometry, self-similarity, only solving for the asymptotic 
structure (which is unconstrained relative to the self-consistent 
global solution that would include the region near the disk), 
treating the accretion disk as a boundary condition, or ignor- 
ing gravity altogether As far as the authors are aware, none 
of the advanced, but still somewhat idealized, relativistic MHD 
jet simulations (Komissarov et al. 2007; Tchekhovskoy et al. 
2009a,b) have been used to study the emission from AGN jets. 
Moreover, in light of the lack of an ab initio understanding of 
particle acceleration within AGN jets, all efforts to model the 
jet emission require ad hoc prescriptions for inserting the emit- 
ting particles. 

Unlike emission maps, the RM maps depend only upon the 
velocity, magnetic field, thermal ion density and electron tem- 
perature. As a consequence MHD simulations are well-suited 
to computing RMs without additional assumptions. Already, 
global axisymmetric MHD simulations have been used to con- 
strain the mass accretion rate and magnetic field geometry of 
the supermassive black hole in the Galactic center (Sharma 
et al. 2007). Computing the RMs of AGN jets requires the 
self-consistent simulating of both the emitting fast jet core and 
the surrounding Faraday screen. However, this is naturally ac- 
complished by GRMHD simulations of jet formation, which 
generally produce a variety of outflow components, including 
an ultra-relativistic jet and a surrounding trans-relativistic wind 
(McKinney 2006). 

Here we present the first analyses of pc-scale RMs associ- 
ated with self-consistent 3D GRMHD simulations of relativis- 
tic jets. For this purpose, we employ the fully 3D GRMHD 
simulation reported by McKinney & Blandford (2009), and 
perform polarized radiative transfer through the 3D volume, 
generating frequency-dependent polarized flux maps and RMs. 
We investigate how these depend upon model, observational 
and physical parameters, paying specific attention to how Fara- 
day rotation within the circum-jet plasma may be distinguished 
from that due to distant Faraday screens. Section 2 summa- 
rizes the computational modeling, including the GRMHD sim- 
ulation, radiative transfer and beam convolution computations. 
Section 3 presents RM maps for a variety of model parameters 
and describes the implications for the location and structure of 
Faraday screen, as well as our theoretical understanding of jet 
formation. Finally, our conclusions are summarized in Section 
4. 

2. COMPUTATIONAL METHODS 

A number of factors affect the computation of the RM maps 
of simulated jets. Foremost among these is the location, struc- 
ture and dynamical state of the Faraday screen itself, which 
we will find to be naturally generated by the circum-jet mate- 
rial immediately surrounding the jet. However, in addition we 
must also carefully account for the relativistic motion within 
the Faraday screen, the suppression of Faraday rotation at 
high electron temperatures, finite radio-beam effects and even 
the location and mechanism responsible for producing the ob- 
served jet emission. In this section we describe how we do 
each. 

2.1. GRMHD Simulations 
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To model the jet system, we use the fully 3D (no as- 
sumed symmetries) general relativistic (GRMHD) simulations 
reported in McKinney & Blandford (2009). We will only sum- 
marize the aspects that are relevant to our present application 
here. The simulations were generated using the conservative, 
shock-capturing GRMHD scheme HARM (Gammie et al. 2003) 
to simulate the evolution of an accretion flow around a rotat- 
ing black hole, and the subsequent self-consistent generation 
of bipolar relativistic outflows. Given an initially dipole mag- 
netic field, embedded within the accretion disk near the black 
hole with spin a/M w 0.92, they have found that generally a 
quasi-stable relativistic jet develops and propagates to large 
radii. While in general many different black hole spins could 
be considered, models with a/M > 0.4 are qualitatively simi- 
lar to those with a/M 0.92 (McKinney 2005). The simula- 
tions run for a duration of about 5 x KfiGM/c^ , by which the 
jet has passed beyond the simulation box, which has a radius 
of lO^GM/c^. The baryon density at the base of the jet, very 
close to the rotating black hole, is determined by limiting the 
comoving electromagnetic energy density per rest-mass energy 
density to no larger than 100, which strictly limits the Lorentz 
factor to less than order 100. The accretion disk has a finite 
extent of about r ^ 40M, which limits the radial range (order 
r ^ lOOM) over which the jet is collimated by the disk, corona, 
and disk-driven wind. This finite collimation extent leads to 
a saturation of the Lorentz factor at roughly 5-10. The effect 
of a more or less extended disk, and so larger or small termi- 
nal Lorentz factors, respectively, can be considered in future 
work. The jet near the black hole is perturbed by disk-driven 
turbulence that induces substructure dominated by the \m\ = 1 
(where m is the azimuthal quantum number) mode at any given 
radius, and there is also substantial power in other spherical 
harmonics. Varying amplitudes of these modes are advected 
with the jet velocity and become part of the jet substructure at 
large radii. A single snapshot at f = 215QGM/c^ was used from 
the simulation data, which is a time corresponding to after the 
jet has fully developed inside the computational box but the 
flow has not significantly reflected off the box. 

The simulations are completely described by the gas rest- 
mass density, p, the gas internal energy, Ug (gas pressure pg = 
(r^- l) Ug), the lab-frame velocities /3, the comoving magnetic 
fields, b, and the spacetime geometry (defined by the metric 
components, g^i,^). While these simulations assume an ideal 
gas equation of state, with a fixed adiabatic index of Tg = 4/3, 
variations in the adiabatic index do not significantly change 
the dynamics of highly magnetized jets (McKinney & Gammie 
2004). From these quantities one can compute the baryon num- 
ber density, n = p/nii, (where nib is the baryon mass) and the 
ideal gas temperature, T = Pg/nk (where k is Boltzmann's con- 
stant). Since we are interested in rapidly accreting, and there- 
fore well-coupled, systems, we assume that the electron and 
ion temperatures are the same. Finally, note that the lab-frame 
and comoving magnetic fields are related by 

B = r(l-/3/3)-b ^ b= ^ (l + r2/3/3) -B, (2) 

where r=(l-/3-) "is the bulk gas Lorentz factor. Thus, for 
large F, e.g., within the jet, B will generally appear toroidally 
dominated for comoving magnetic fields with even moderate 
pitch-angles. 

While the black hole metric produces natural length 
(GM/c^) and time (GM/c^) scales, the GRMHD simulations 
have no natural scale for any density, such as the plasma den- 



sity, gas pressure, and magnetic field strength. Thus we are 
free to renormalize the simulation data by choosing a single 
arbitrary number, which we take to be the black hole mass ac- 
cretion rate, determined at the horizon by, 

M= f dn^pu', (3) 
Jll, 

where g = Det(g^y) is the determinant of the metric, u'' is 
the contravariant radial 4-velocity and r/, is the radius of the 
horizon. This gives a mass scale of MGM/c^, or equiva- 
lently, a density scale of Mc^ /G'M^. In the following sec- 
tions, this will be chosen to qualitatively reproduce observa- 
tions, and will vary between 10~^-10~^ times the Eddington 
rate (defined in terms of the Eddington luminosity, Lndd, to be 

MEdd ^ lOLEdd/c"). 

2.1.1. Extrapolation 

Characteristic black hole masses for the AGNs of interest 
range from IO^Mq to IO'^Mq. Given that we are seeking to ad- 
dress the Faraday rotation measure distributions of AGN jets on 
pc scales, we are concerned with the jet structure on scales of 
10*'GM/c--10'*GM/c^, respectively. Even at the highest mass 
range, the simulated region is too small. Propagating the jet 
to a distance of lO^GM/c^ within the GRMHD simulation is 
presently computationally infeasible. Thus we are either faced 
with employing unphysically large masses (~ 1O"M0) or ex- 
trapolating to larger distances. In this paper we will do both. 
Here we describe the procedure we employ for the latter, justi- 
fied by the asymptotic, power-law behavior of relativistic out- 
flows. 

Beyond the fast magnetosonic surface, highly magne- 
tized, relativistic jets reach the so-called monopole phase 
(Tchekhovskoy et al. 2009a,b). Once the jet is in this phase the 
azimuthally and time-averaged fluid quantities are very nearly 
self-similar and the jet moves along conical flow lines, follow- 
ing power-laws for the rest-mass density, internal energy den- 
sity, components of the velocity and magnetic field, and a log- 
arithm for the Lorentz factor (McKinney 2006; Tchekhovskoy 
et al. 2009a). As noted in studies of otherwise similar 2D jet 
simulations that extend to larger sc ales (lO^GM/c") than the 
present 3D simulations, the location of the fast magnetosonic 
surface is at r ~ lO^GM/c^ (McKinney 2006), suggesting that 
we may approximate the structure of the outflow at large radii 
by advecting the simulated quantities outside of the fast mange- 
tosonic surface. This is borne out by inspection of the time and 
azimuthally averaged jet quantities in the full 3D simulation as 
well'. Extrapolating the time-dependent, 3D fluid quantities 
assumes that the spectrum of fluctuations is also self-similar, 
though this is natural in the context of turbulence. 

In axisymmetry, the self-similar structure of the fluid quan- 
tities is constrained by a number of MHD invariants. These 
are described in detail in McKinney (2006), and therefore only 
summarized for the conical regime here. As a consequence 
of these constraints, we must define only one free parame- 
ter, corresponding to the rate of dissipation, which is naturally 
bounded. Generally, within the conical regime, the fluid quan- 
tities are advected radially, and are thus related to their values 
along radial flow lines. In particular, the radial dependencies 

' While the 3D GRMHD simulation does appear to follow the same generic 
stracture found in the earlier large-scale 2D simulations reported in McKinney 
(2005), it has insufficient radial dynamic range to unambiguously define the 
relevant asymptotic behaviors of the fluid quantities on its own. 
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(at fixed 0) of the toroidal and poloidal lab-frame magnetic 
field components are roughly given by B'^ oc (rsin6') 



1 -1 
oc r 



and B'' (xB'f'/ 



respectively . Within the coast- 



ing regime the Lorentz factor grows with radius only as a 
small power of a logarithm (i.e., roug hly (ln(r))i/3), and thus 
we approximate it as fixed, giving c>c 1 along radial lines. 
The toroidal velocity is constrained via angular momentum 
conservation, giving /^"^ cx (rsin6')~' cx r"' along radial lines. 
Mass conservation then implies that n cx B^ /u'' cx r~- along 
radial lines. Defining the structure of the internal energy de- 
pends critically upon the efficiency with which magnetic dis- 
sipation heats the baryons and electrons. In principle, this 
is measured in the simulations, with McKinney (2006) find- 
ing u cx r~'^ along flow lines. However, this cannot continue 
indefinitely; the magnetic energy density available decreases 

as (where b~ = |b|~ - (b • /3)^ is the square of the mag- 

netic field strength in the comoving frame, not to be confused 
with simply |b|^), placing an upper limit upon the magnetically 
heated baryon energy density. A lower limit can be found by 
assuming that the flow is adiabatic, giving u cx ii^/^ cx r~^l^ . 
We consider both cases, adiabatic expansion of the flow as 
well as setting u cx r"^^ until it reaches e\r jZ-K, for e ~ 0.01- 
1 .0. In summary, within the self-similar, conical regime, the 
power-law indices of «, /3 and B are determined explicitly. The 
specifics of magnetic dissipation can effect the internal energy 
of the baryons, and thus we consider two limiting cases: adia- 
batic and fully magnetically-thermalized flows. 

In practice, to perform the extrapolation we define a 
piecewise-linear remapping. At locations within Ri^^ = 
IQ~ GM/c~, roughly the location of the fast magnetosonic sur- 
face, we do no extrapolation, returning the simulation data as 
is. This ensures that we do not extrapolate the non-asymptotic 
jet structure to large distances as well as restricting the extrap- 
olation to the outflow (since the accretion disk is located well 
within this radius for the duration of the simulation). Outside 
of /?fnis we linearly stretch the coordinates so that the outer sim- 
ulation boundary, R^im, is now located at R^ax- That is, the rela- 
tionship between the radius in the physical domain, r, and that 
within the simulation, x, is given by 



if r < 7?fms 

r» — (r-Rfms) otherwise. 



In terms of this, following the discussion above, we set the 



McKinney (2006) found that the jet can shghtly de-colhmate at large radii 
such that, e.g., oc r"'^ for some field lines. However, much of the jet 
power flows along those field lines that have B*^ oc r~K Also, we consider 
a single simulation with ci/M fa 0.92, while models with much smaller a/M 
produce a black hole jet that only starts to have a doininant lab-frame toroidal 
field beyond the Alfven radius (hght cyhnder for pulsars) located at roughly 
4cr;,/(a/M), which can be quite large for small a/M. However, below a/M 
0.4 black hole diiven jets are dominated by the Keplerian disk-driven wind 
(McKinney 2005). As a consequence, for low black hole spins the disk driven 
wind overpowers any observations of the black hole driven jet. Therefore, the 
single a/M ^ 0.92 simulation exemplifies any model that would focus on black 
hole driven jets. 



values of the fluid quantities in terms of the simulation values. 



B. 



{r,e,cj,)=[^BHx,0,^) 
BL(r,0,0)=(^)'BV,e,0) 



(5) 



/3f,,(r,0,0) = /3^(x,0,0) 
«ext('", 0,4}) =(^^^ n{x, 9, 0) 

/ X\ 

and subsequently define b given Equation (2). We consider 
fl« = -8/3 and -1.3 (limiting it from above by elp- jZ-K), corre- 
sponding to the lower and upper limits upon the magnetic dis- 
sipation rates, respectively. Henceforth we will drop the sub- 
scripts upon the fluid quantities, though it should be understood 
that in all cases we are referring to the extrapolated values. 

From Equation (5) it is clear that there are important differ- 
ences between the pc-scale structure of low-mass (1O**M0) and 
high-mass (10'"Mq) AGNs. Generally, we expect that at the 
same physical distance the magnetic fields of low-mass AGNs 
will be even more toroidally dominated (in all locations) than 
their high-mass analogs. For similar reasons, the velocities of 
low-mass AGNs will be more radially dominated than those of 
high-mass AGNs. 

Our results are primarily dependent upon the conically ex- 
panding region beyond r ^ lOOGM/c", except when the black 
hole mass is chosen to be unphysically large. For realistic 
black hole masses and pc-scales, the jet magnetic field becomes 
toroidally dominated in both the lab and comoving frames. 
Jets that are toroidally dominated within the comoving frame 
are commonly claimed to be unstable to helical kink (screw) 
modes (see, e.g., the discussion in Lyutikov et al. 2005 regard- 
ing toroidally dominated cylindrical jets). However, the conical 
jets considered here are stabilized by lateral expansion, result- 
ing in causal disconnection across the jet, quenching the growth 
of unstable modes (McKinney & Blandford 2009). Thus, the 
substructure of conical jets at large radii should be indicative 
of the interaction between the black hole, accretion disk and jet 
close to the jet launching point, which has subsequently been 
advected to large radii without significant additional disruption. 
This is consistent with the above extrapolation procedure that 
essentially advects the jet substructure to large radii. However, 
we caution that such a re-scaling cannot be universally applied 
to AGN jets, failing for those jets that may be significantly af- 
fected by their environment (e.g., those associated with FRTs). 
In addition, while this implies the jet is roughly stable, some 
flaring behavior due to small-scale MHD instabilities may also 
be expected (Sikora et al. 2005). 

2.2. Rotation Measures 

Provided the thermal electron temperature, proper density, 
four- velocity and magnetic field strength and geometry, we can 
now compute the RM along any given line of sight. Following 
Broderick & Loeb (2009), we do this directly via 



RM = 0.812 X 10' 



k-/3) -bc/fradm"-, (6) 
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Figure 1. Schematic representation of the relevant portions of the simulated 
jet- wind-disk system and their relationship to example lines of sight. The ther- 
mal particle proper density is shown by the color map, extending over 9 orders 
of magnitude. The disk and jet regions are delineated by the thin and thick 
white lines, respectively. In between is the corona and wind (including both 
bound and unbound components). The limited extent of the disk and wind is 
due both to finite initial extent of the disk and to rotational collimation. The 
observer is shown in the lower-right, with the image plane immediately above. 
Lines of sight along which the RM and flux are integrated are shown, with the 
former being integrated along only the red portion. Dotted lines of sight pass 
through the initial transient portion of the simulation and are thus discarded 
(see Section 2.5). Finally, the definition of is presented explicitly. 

where g = r(l-/3 k) and we have supplemented this with a 
finite temperature term^, given by, 

1 7-1 / kT \ 
f{T)= -r + ——log{-/) where 7 = max 1, r , (7) 

7^ \ nieC^ J 

where 7 is the thermal particle Lorentz factor, not to be con- 
fused with the bulk Lorentz factor of the out-flowing gas 
(Huang et al. 2008; Shcherbakov 2008). Note that the RM, 
as we have defined it, depends solely upon quantities that have 
been simulated, with the integration geometry shown in Figure 
1 . All integrals are performed in flat space, ignoring any con- 
tributions from gravitational lensing, an approximation that is 
well-justified at the distances responsible for the resolved po- 
larized emission. With this restriction. Equation (6) is correct 
for relativistically moving Faraday screens. 

The integral in Equation (6) is performed from a vertical 
mid-plane, to a distant image plane, orthogonal to the line of 
sight. In principle, this overestimates the observed RM because 
it ignores depolarization along the line of sight resulting from 
Faraday rotation within the emitting region. However, in prac- 
tice the emitting and rotating regions are well-separated, as il- 
lustrated explicitly in Figure 2, which shows the fractional local 
contributions to the Faraday rotation and emission (see the fol- 
lowing section) along a representative line of sight. That this 

^ The rate of Faraday rotation depends upon the nature of the underlying 
electromagnetic eigenmodes of the plasma. For cold plasmas these are cir- 
cular, reproducing the standard result. For ultra-relativistic plasmas these are 
substantially elliptical, strongly suppressing any Faraday rotation. 




/ (GM/c2) 



Figure 2. Fractional integrated intensity (black) and rotation measure (red & 
blue dashed lines, corresponding to positive and negative RMs) as a function 
of distance as measured from the jet mid-plane along a number of exemplar 
lines of sight, which were chosen such that they span the width of the jet. After 
reaching 99% of its final value, the intensity is shown by a dotted line. In 
all cases, along any given line of sight the region responsible for the RM is 
well-separated from that responsible for producing the observed emission. 

is the case is also supported directly by the observation of sig- 
nificant jet polarization fractions, ruling out substantial Fara- 
day rotation within the emitting region (Burn 1966; Gardner & 
Whiteoak 1966). To generate an RM map, we then repeat this 
process many times at different positions on a rectilinear grid 
located on the image plane. Examples of the resulting maps 
can be found in Figure 3 and Section 3. 

2.3. Jet Synchrotron Emission 

While we may now compute the RM along any line of sight, 
given a complete description of the circum-jet material, in prac- 
tice RMs can be measured only for those lines of sight along 
which the jet is sufficiently bright. As a consequence, in order 
to compare with observations of AGN jets, we must also model 
the jet emission. This will serve two purposes: first as a mask 
upon the RM-maps, identifying the regions that are observable, 
and second as a weight for determining the beam-convolved 
RM maps, as described in Section 2.4. 

2.3.1. Particle Acceleration 

Arising from nonthermal synchrotron emission within the 
jet itself, the emission is much more poorly constrained by 
the simulations. The difficult arises primarily from the lack 
of an ab initio understanding of the mechanism by which the 
radiating nonthermal electrons are accelerated. Here we con- 
sider a simple, qualitative model, appropriate when the ac- 
celerated electrons are the result of the dissipation of MHD 
turbulence and/or magnetic reconnection within the highly- 
magnetized jet"*. In this case the reservoir of energy that is 

Internal shocks have also been suggested as a mechanism for producing 
the observed nonthermal particles. We have considered a variety of shock- 
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Figure 3. Ideal-resolution and beam-convolved polarized flux with unrotated polarization ticks (upper-left), spectral index (upper-right) and RM maps (lower-left) 
for our canonical jet model as observed at 15 GHz (see Table 1 for a complete list of relevant parameters). In addition the contributions from the two terms in 
Equation (21) are shown expHcitly (lower-right). In all panels total flux contours are overlaid, in logarithmic factors of 2, with the minimum contour corresponding to 
0.98mJymas"^, and all other quantities are clipped at a polarized flux of 1 mJymas"-. Beam-convolved maps show the beam size (0.3 mas FWHM) in their lower-left 
corner. 



available to be tapped via particle acceleration is that associ- 
ated with the electromagnetic field. Thus, we set the comoving 
energy density in nonthermal electrons, Wnth, proportional to 
that of the simulation electromagnetic field, b^/diir. This as- 
sumption is consistent with minimum energy arguments used 
to infer the expected density of emitting particles and mag- 
netic field strengths for emission due to synchrotron emission 
within AGNs (Burbidge 1956). We further restrict the parti- 
cle acceleration by requiring that it occur only in magnetically 
dominated regions, i.e., where b~ /%tt > pc^, which ensures that 
particles can be efficiently accelerated to large Lorentz factors 

acceleration schemes, based upon identifying the locations of resolved shocks 
within the simulation and setting the nonthermal energy density to a shock- 
strength dependent fraction of the post-shock internal energy. These generally 
result in disk-dominated emission, and are incapable of reproducing the ob- 
served extended jet emission. This implies that if shocks are responsible for 
producing a substantial portion of the nonthermal electrons, they must be mi- 
croscopic or produced by an interaction between the jet and a suiTounding 
ambient medium that is not included in the simulation we consider 



via MHD turbulence and/or magnetic reconnection (Lyubarsky 
2005; Lyutikov 2006). That is, explicitly, we set 



"nth = r] 



47r 



1 -l-exp 



pc- 



(8) 



and 77 is the efficiency with which magnetic energy is con- 
verted into that in nonthermal particles [recall that = |b|^- 
(b /3)^]. In magnetically dominated regions this reduces to 
"nth = rib^/^TT, and is exponentially suppressed otherwise. This 
is motivated first by the physical assumption that conver- 
sion of electromagnetic energy into nonthermal particles is 
less efficient when the electromagnetic energy per particle is 
low (e.g., relativistic magnetic reconnection requires, roughly, 
b^/{8n) > pc^, Lyutikov & Uzdensky 2003). However, in 
practice relaxing this suppression results is broad flux distri- 
butions that are inconsistent with observations of AGN jets. 
As a physically motivated alternative, one might choose Mmh 
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to be proportional to the square of the comoving current den- 
sity, indicative of where current sheets are dissipating (Lyu- 
tikov et al. 2005). However, the current density in the simu- 
lations is unlikely to be properly resolved and outside the cur- 
rent sheet the properties of the current layer are determined by 
the comoving electromagnetic energy density, suggesting that 
Mnth C!C also roughly tracks the importance of current sheets 
if they were present. We leave the consideration of other pre- 
scriptions for the nonthermal particle density for future work. 

We further assume a power-law electron distribution, with 
index -2 (resulting in a spectral index of 0.5) and minimum 
particle Lorentz factor, jmin, of 10". Since our primary con- 
cern is to identify bright and dim regions, moderate changes in 
either of these do not significantly change our results. 

2.3.2. Polarized Radiative Transfer in Bulk Flows 

Due to the strong magnetic field within the jet and the rela- 
tively low density, the jet emission will be dominated by syn- 
chrotron. On pc scales, AGN jets are generally optically thin, 
while the bright radio core is optically thick. Thus, synchrotron 
self-absorption is critical to obtaining total jet fluxes and mod- 
elling the jet spectral index distribution. Furthermore, in order 
to compute the beam-convolved RMs we will need intrinsic 
Stokes /, Qo and Uo- Thus we are necessarily faced with com- 
puting the self-absorbed, polarized emission. 

The transfer of polarized radiation through relativistic envi- 
ronments has been discussed by a variety of authors and thus 
here we only summarize the procedure (see, e.g., Lindquist 
1966; Connors & Stark 1977; Connors et al. 1980; Laor et al. 
1990; Jaroszynski & Kurpiewski 1997; Bromley et al. 2001; 
Broderick & Blandford 2004; Beckwith & Done 2005; Brod- 
erick 2006; Broderick & Loeb 2006). This is simplified some- 
what by the assumption that within the regions of interest we 
may approximate the spacetime as being flat, requiring only 
that we consider the relativistic motion of the emission region. 
This results in the Doppler boosting and beaming of the emis- 
sion as well as relativistic aberration of the polarization. All 
of these can be taken into account by noting that the photon 
occupancy numbers, cx {I,Qq,Uq)/v^, are Lorentz scalars, al- 
lowing us to rewrite the standard radiative transfer equation in 
the emitting frame. From this we obtain. 



leaves the spectral break associated with 7min, giving 



IN 



1 Cg Cv 



Go|=/. I eel-a. I Ce 1 one 
uJ VuJ \Cu 1/ 



(9) 



where j^, = jgv/g~ and oii, = gagi, are Lorentz covariant 
analogs of the lab-frame emission and absorption coefficients 
{ji, and a^, respectively), which with equ and Cq.u, are dis- 
cussed in more detail below. 

Generally, the form of j,y depends upon the locations of the 
minimum and maximum nonthermal particle Lorentz factors 
(we have specified only the former), and the location of a cool- 
ing break (located where particle acceleration balances syn- 
chrotron cooling). However, the cooling time at radio wave- 
lengths exceeds the jet propagation time at pc scales, and thus 
we may safely ignore cooling in the computation of the pc- 
scale radio images. Furthermore, while the maximum particle 
Lorentz factor generally depends upon the details of the par- 
ticle acceleration mechanism, typical estimates, given the jet 
magnetic field strength, are much larger than the Lorentz fac- 
tors of the particles responsible for the radio emission. This 



- „ "nth . 1 5 Q 

J^ = Cj-^ UBSm -d 



gVmin 

0.5 



0.5 



1/3 



otherwise. 



(10) 



in which Cj ~ 60Jycm^erg 'Hz ' sr vg = eb/lirmgC is 
the cyclotron frequency in the comoving frame, sint? = 



g^b^- [{k-0)-h\' I gb is the sine of the angle between the 

line of sight and the magnetic field in the comoving frame, and 
t'min = 7^;jjfBsin'(5/g is the frequency associated with spectral 
break due to 7niin. Given j^, we may determine via Kirch- 
hoff 's law, giving. 



Mnth sin" I? 



7min Vb 



/ \ 3 / \ 1/3 

\gVn 



mm / \ ^mm 



KVrr 



otherwise, 



(11) 



where Ca ~ 1.3 x lO^'cm"'. 

The polarization terms are affected by both the relativistic 
aberration associated with the bulk motion of the jet and the 
orientation of the magnetic field. Both of these can be ad- 
dressed by parallely propagating an orthonormal tetrad along 
the ray trajectory, defined at the image plane such that two of 
the basis vectors are aligned with the directions defining the 
observed Stokes parameters (and thus orthogonal to the line of 
sight). In principle, this basis may then be boosted into the 
emitting frame, directly accounting for relativistic aberration, 
and then compared with the projected direction of the magnetic 
field. In practice, it is sufficient to define the angle, ^, between 
the initially vertical basis vector (associated with Q) and the 
projected magnetic field direction, in the emitting frame. As 
shown in Appendix A, this is defined up to a sign by 



cos^ = 



z- (1-kk) -(b + r/S/jcosi?) 



(12) 



1 - (z • k) Zjsini? 



The sign of f may then be obtained by a similar construction 
obtaining sin^ from the relationship between the initially hori- 
zontal basis vector and the projected magnetic field in the emit- 
ting frame. However, it is sufficient to note the direction of b in 
the lab frame since boosting preserves the ordering of vectors, 
e.g., if b is to the right of z (as seen by looking down the line 
of sight) in the emitting frame, then it will be so in all frames. 

Thus, sin^ = sgn [b • (z x k)] \J 1 -cos^^. In terms of these, 
the polarization terms are given by 



9 9 

eG = -Y^cos2^, ej/=-— sm2^, 

3 3 

Cg=-4Cos2^, Cf/=-4sm2^, 



(13) 



where the particular numerical coefficients are set by the choice 
of the electron spectrum index and the overall sign reflects the 
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fact that in the emitting frame the polarization is orthogonal to 
the projected magnetic field. 

As with the expression for the RM, Equations (9-13) are 
Lorentz covariant. Equation (9) is integrated throughout the 
emitting region along a given line of sight ending upon the im- 
age plane. However in practice, the emission is limited to the 
jet and associated accretion disk. Again, this procedure is re- 
peated at a number of points upon the image plane, producing 
an image of the jet itself. Upon choosing a black hole mass, M, 
distance and beam size, we generate a normalized radio map of 
the polarized jet emission, examples of which may be found in 
Figure 3 and Section 3. 

In the following section we will find that the spatial distri- 
bution of the spectral index of the polarized emission, «£, will 
be important in determining the observed RM maps. As a con- 



sequence, we also construct a/, = -d log y QI + Uq /d log ly and 

a = -dlogl /dlogi/ by differencing polarized emission maps 
at nearby frequencies. Both ai and a (an example of which 
is shown in Figure 3) are generally core dominated with a size 
indicative of the accretion rate. 

2.4. Beam Convolution 

Radio maps of AGN jets are generally under-resolved, with 
the consequence that a great deal of the jet substructure is 
smoothed out by the radio beam. However, in the case of the 
RM maps, we must take some care since RM is not observed 
directly. In this section we describe how we perform an appro- 
priate beam convolution operation on the RM maps. We be- 
gin with a general discussion regarding finite beamwidth RM 
maps, and then move onto the particulars of the scheme we 
utilize in this paper 

2.4.1. Finite Beamwidtli RM Maps in General 

It is possible to infer the beam-convolved RM maps without 
producing fully simulated observations, i.e., without explicitly 
producing polarized images of the jet at multiple wavelengths. 
This is done by careful considering what is, and is not, directly 
measured. In practice, the RM is a derived quantity, computed 
from directly measured Stokes parameters Q and U. Along 
with /, these are limited by the effective telescope beam, with 
the observed map being related to the infinite-resolution values 
via a convolution: 

{i)in) = J dnBin'-n)iin'), (i4) 

with similar relations for Q, U and V, in which B{fl) is the 
beam profile (which we will ultimately take to be Gaussian). In 
contrast, since RM is non-linearly related to Q and U, the rela- 
tionship between the finite-beam and the ideal-resolution RMs 
is not a simple convolution. Nevertheless, as we show below, it 
may be written as the combination of weighted convolutions. 

The effect of an intervening Faraday screen is to rotate Q into 
U and vice versa: 



2 = Lcos (2*o + 2RMA2 
U = Lsm (2*o + 2RMA^' 



(15) 



where L= y Q^ + U^ is the intrinsic (ideal-resolution) polarized 

flux, 2o =i'Cos(2\E'o) and f/o = Lsin(24'o). The observed Q and 
U are then 

(e) = (Lcos(2*o + 2RMA2)) 
(f/) = (Lsin(2*o + 2RMA-)) , 



and is now explicitly dependent upon the spatial distribution of 
the RM and A. The observed polarization angle is 



*obs = 4 arctan 



(U) 



2 V (2) 

This implies that the measured RM is then 



(17) 



d\^ 



1 



dX^ 



Ignoring the wavelength dependencies in the beam size 
(which are explicitly accounted for with match-resolution im- 
ages, produced by convolving high-frequency images with the 
lowest frequency-interval beam) d/dX- commutes with the 
beam convolution. Thus, with Equations (15), 



(19) 



where we have assumed that ^Pq is wave-length independent 
(i.e., Faraday rotation is the only frequency dependence in 'i/). 
Similarly, 



(20) 



As a consequence. 



RMnhs = 



{Q){QRM) + {U) (URM) 

{Qf+{uf 

, {Q){Uc^L)-{U){QaL) 
4X^({Qf+{Uy 



(21) 



The first term corresponds to a weighted average of the RM 
across the beam. Unresolved correlations between the polar- 
ized flux and the RM can produce substantial discrepancies be- 
tween the intrinsic and observed RMs. At this point we can 
identify some of the important differences between observa- 
tions of RMs and the polarized fluxes. If the RM is fixed across 
the beam we have RMobs = RM trivially. This is true even if 
there are rapid variations in the intrinsic polarization^. The 
converse is not true; order Itt / X^ variations in RM across the 
beam results in the well-known Faraday beam depolarization. 
More important are correlated variations in the polarized in- 
tensity and the RM. Since RMobs is a weighted average, clear 
large-scale structures in the RM can be obscured by the pres- 
ence of beam-scale structure in L. Thus even a uniform RM 
gradient can appear non-uniform at the beam-scale when illu- 
minated by a spatially variable source, and in particular near 
the radio core where all quantities are changing on scales small 
in comparison to the observing beam. 

The second term is associated with spatial variations in the 
polarized spectral index, vanishing identically when is con- 
stant. This is due to the potential for frequency dependent 
competition among unresolved regions with different intrin- 
sic polarizations. Spatially varying spectral indices have been 
found in a number of sources. Most common are swings from 
ai ~ 1-0.5 within the optically thick AGN core and settling to 

^ Though in this case the measured polarized flux may be so small as to 
make measuring the RMobs no longer feasible. 
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ai ~ -0.5 at angular distances of more than ~ 1 mas, where 
the jet is optically thin, though cases where similar magnitude 
variations across the jet occur are known, presumably due to in- 
teractions with the AGN environment (O' Sullivan & Gabuzda 
2009b, a). For comparison with these we also compute the 
beam-convolved spectral index; 



_ d\og{I) ^ {la) 
dlogX (I) 



(22) 



Examples may be found in Figure 3 and Section 3. 

2.4.2. Beam-Convolved RM Gradients 

We will remark upon the implications of finite beams for the 
RM-maps of AGN jets in detail in the following sections. How- 
ever, a simple, qualitative example directly illustrates the im- 
portance of beam convolution upon the search for, and study 
of, transverse RM gradients within AGN jets. We will find that 
generally it is crucial to resolve the transverse structure in both 
the emission and the RM distribution in order to obtain accu- 
rate measurements of the true RMs. 

We briefly consider a uniformly polarized ID, cylindrically 
symmetric jet with a Gaussian transverse flux profile with 
width (jj, a circular Gaussian beam with width cjb, and a lin- 
ear RM-profile. The typical magnitude of the RMs under con- 
sideration is lO^radm"^, which at 15 GHz produces Faraday 
rotations of ^ 0.4rad. For the present purpose this is suffi- 
ciently small that it may be safely neglected (though we have 
confirmed this via explicit computation). Thus, without loss of 
generality we may set. 



RM = RM^);c, B = 



Ittctb 



Q = 



(23) 



where we have rotated the axes defining the Stokes parameters 
to force U to vanish. With these the various beam-convolved 
quantities may be computed analytically, giving 



dyB{y-x)Q(y)-- 



(e>w= 

(eRM)(.^)= / dyB{y-x)Q{y)mA{y) = 

J—OG 

which then implies 



-X- /2{a]+ul) 



RM, 



obs 



:RM 



(24) 



(25) 



That is, for linear RM gradients, the beam-convolved RMs are 
generally smaller than the true values, and is a relatively strong 
function of the ratio of jet width to beam size. As a conse- 
quence, when the transverse jet structure is not resolved, even 
apparently linear RM gradients will not be accurate measures 
of the true jet RMs. We will find that Equation (25) will be 
useful in interpreting the significance of unresolved, but never- 
theless robust, RM gradients. 

2.4.3. Computing Beam-Convolved RM Maps in Practice 

Where the gradients in RM are sufficiently high, the proce- 
dure by which we have computed the RM, {I,Qo,Uo) and ai 
maps fails to resolve the rapidly varying polarizations. That 
is, when the image grid scales (generally much smaller than 
the beam size) are Afl^ and Afly, if Ail^dRM/dil^ ^ or 



A^lydRM/dfly > TT, we must supplement the inferred Q and 
U maps with a sub-grid model. Since the grid scale is much 
smaller than the typical beam scale, we may assume that 8(0,) 
is roughly constant over a single pixel, and thus the primary 
effect corresponds to a grid-scale analog of beam depolariza- 
tion. In this case, it is straightforward to show that the effective 
polarization fraction is 



. - ^ t/RM ,\ . / ^ dRM^, 

mi = mn) sine Aitjc A smc Ally A" 

dilx ) \ ' dily 



(26) 



Thus, the procedure for producing the a finite-beam RM map 
is to first compute the ideal resolution RM, {I,Qo,Uo) and ai 
maps, pixel-depolarize these maps using Equation (26), and 
perform the convolutions in Equation (21) assuming a Gaus- 
sian beam. Examples of the resulting RM maps are presented 
in Figure 3 and Section 3. 

2.5. Geometric Limitations 

There are necessarily geometric limitations upon the jet in- 
clinations, 0, for which we can produce RM maps. For <d less 
than the jet half-opening angle, ~ 3° at Rmax, we are looking 
down the jet directly. As a consequence, the resulting RMs are 
indicative of initial transients in the simulation, and of the radio 
lobe in real systems. Thus, we restrict ourselves to > 10°. 

However, these transients combined with the finite simula- 
tion time constrain the angular size of the image region more 
generally. The projected length of the jet is roughly Lpioj ~ 
^maxSinB. Lines of sight which pass through regions more 
distant than this, shown by the dotted lines in Figure 1, are nec- 
essarily indicative of the initial transient period. Thus, we limit 
ourselves to lines of sight that pass through the non-transient 
portions of the jet, and therefore restrict the region on the image 
plane we ultimately consider. Since this constraint becomes 
more severe at small Q, becoming impossible to satisfy when 
8 is equal to the half-opening angle of the jet, the credible 
portions of the RM maps are necessarily different sizes for dif- 
ferent inclinations. 

3. JET RM MAPS AND THEIR IMPLICATIONS 

In this section we present RM maps for a variety of model, 
observational and theoretical parameters. These include varia- 
tions in AGN model parameters (e.g., the jet orientation, black 
hole mass and magnetic dissipation rate; §3.1), observing pa- 
rameters (observation frequency and beam size; §3.2) and the 
structure of the Faraday screen (§3.3). To ensure that we are 
comparing physically similar systems, we define a canonical 
model: A central black hole with mass 3 x 1O**M0 accreting 
at 1% of the Eddington rate located at a distance of 300 Mpc 
and viewed at an inclination of <d = 20°. For the extrapola- 
tion we choose Rnnx = ^O^Rsim = lO^GM/c^ ~ 15pc, and select 
flu = -8/3, corresponding to the adiabatic advection model. Fi- 
nally we choose 77 = 0.02 such that the jet flux is roughly 4 Jy 
at 15 GHz. These choices are not unique, and we will consider 
explicit departures from this canonical model in all respects; 
a list of the parameters of the models we present are listed in 
Table 1. 

Figure 3 illustrates the general morphologies of the flux, po- 
larization, spectral index and RM associated with our simu- 
lated jets. As in all subsequent RM maps, we show the polar- 
ization properties only where they are likely to be observed, 
requiring the polarized flux exceed 1 mJymas"^ ^. Specifically, 

^ We choose a limit on a fixed surface brightness, as opposed to Jy per beam. 
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Table 1 

Jet model parameters 



Figure v beam" M Mr] a,, e ^isgHz 

(GHz) (mas) (R,;,,) (lO^Mo) (Mfidd) (Jy) 

3,4,6.7,10,12 15 0.3 20° 10^ 3 0.01 0.02 -8/3 - 3.83 

532 15 03 10° H? 3 0.005 002 ^8/3 - 3.77 

5,12 15 0.3 30° 103 3 0.018 0.02 -8/3 - 3.89 

6 15 03 20° 10^ 30 002 ^8/3 - 4.00 
6,12 15 0.3 20° 10 300 10"* 0.02 -8/3 - 6.67 

6 15 0.3 20° 1 3000 10'^ 0.02 -8/3 - 12.7 

702 15 03 20° W 3 OOl O02 ^3 1 3.83 

7 15 03 20° 10^ 3 0.01 0.02 -1.3 0.1 3.83 

7 15 0.3 20° 10^ 3 0.01 O02 -1.3 OOl 3.83 

9 5 03 20° H? 3 OOl O02 ^8/3 - 3.35 

9 8 0.3 20° 10^ 3 0.01 0.02 -8/3 - 3.57 

9 22 0.3 20° 103 3 0.01 0.02 -8/3 - 3.97 

9,12 43 0.3 20° 10^ 3 0.01 O02 -8/3 - 4.15 

8 15 O04 20° H? 3 OOl O02 ^8/3 - 3.83 
8 15 Ol 20° 103 3 0.01 O02 -8/3 - 3.83 

8,12 15 09 20° 10^ 3 0.01 O02 -8/3 - 3.83 



"FWHM 

in all cases we see an optically thick, small core that domi- 
nates the total and polarized emission (see the left-most and 
right-center upper panels). The ordered magnetic field within 
the jet results in an ordered polarization orientation, departing 
from the preferred direction only near the jet edge due to aber- 
ration and projection effects. Similarly, the ordered magnetic 
field in the material surrounding the jet results in an strongly- 
ordered RM-map, exhibiting nearly linear transverse gradients 
along the length of the jet, and peaking near the core (lower-left 
panel). 

These features, however, are somewhat obscured upon con- 
volving with an observing beam. While generally the optically 
thick core is unresolved, it dominates the spectral index within 
a beamwidth of the core. The large RMs near the core result 
in substantial beam depolarization, shifting the peak polarized 
flux away from the peak total flux. While the derotated po- 
larization vectors retain their order, the polarization map can 
be dramatically effected, exhibiting the nearly linear gradients 
only when located more than a beamwidth from the core. This 
is primarily due to the first term in Equation (21), with the sec- 
ond term typically making ~ 10% corrections. While this sug- 
gests that the spectral index terms may be ignored, it should 
be noted that they produce gradients in random directions, and 
therefore directly obscure the transverse gradients indicative of 
the ordered toroidal magnetic fields. 

3.1. Dependence of the RM upon the AGN Parameters 

Here we address the importance of the jet and black hole 
parameters upon the RM maps. These arise via a number of the 
facets of our jet model, depending upon both the fundamental 
properties of the AGN itself as well as the assumptions that 
enter into the extrapolation scheme. 

3.1.1. Jet Orientation 

Chief among the jet properties affecting both the RM maps 
and the emission is the orientation. Having a self-consistent 
3D-jet model provides the unique opportunity to assess the 
systematic uncertainty associated with non-axisymmetric jet 
structures, providing a fundamental lower limit upon the un- 

so that we may later compare maps with dilferent beam sizes on equal footing. 



certainty associated with efforts to infer information about the 
circum-jet material from RM distributions. 

Figure 4 shows the RM maps for our canonical jet model 
as viewed from a variety of azimuthal angles, $ (though all 
at the identical inclinations). To the right of each RM map 
we show example transverse profiles, taken near the core, far 
from the core and at an intermediate location. Within these the 
transverse profiles for views with $ = 0°, 45°, 90°, 135°, 180°, 
225°, 270° and 315° are shown in addition to that associated 
with the neighboring map, providing a sense of the variations 
along each transverse section. Large variations among views 
indicates large inherent uncertainties in the RM values and their 
gradients due to the 3D jet structure. 

Qualitatively, all azimuthal views are similar, exhibiting sub- 
stantial beam depolarization in the radio core, and a region of 
large RMs near the polarized flux maximum. Despite being 
nearly linear at all locations shown prior to beam convolution, 
where the transverse jet structure is not resolved the RM pro- 
files vary significantly. However, despite the large structural 
variations in the jet itself, at positions further than a single 
beam, the transverse gradients become much more robust and 
are indicative of those present at ideal resolution. 

Note that while the presence of substantial variations be- 
tween transverse RM profiles does directly imply large uncer- 
tainties in the RM gradients, the inverse is not true. That is, 
consistent RM gradients do not necessarily imply that these are 
indicative of the circum-jet structure. As implied by Equation 
(25) and as we shall see explicitly below, it is possible for large 
beams to artificially produce essentially linear gradients that 
are considerably smaller than the associated intrinsic values. 

Jet inclination is known to have dramatic consequences for 
the appearance of AGN jets. This is also true for jet RM-maps. 
Generally, where the emission becomes increasingly core dom- 
inated the finite beam effects (see Section 3.2) are exacerbated. 
In addition, at different jet inclinations the RMs probe differ- 
ent circum-jet regions, with smaller O being sensitive to larger 
radii. 

In Figure 5 we show the jets with inclinations of = 10° and 
30°. In order to ensure we are comparing observationally sim- 
ilar systems, we adjust the accretion rate so that the total flux 
is comparable to our canonical model, requiring variations of 
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Figure 4. RM maps at 15 GHz for our canonical model (see Table 1 for a complete list of relevant parameters), as seen from a number of azimuthal angles, <E> 
(listed in the lower-right of each RM map). For reference total flux contours are overlayed, in logarithmic factors of 2, with the minimum contour con'esponding to 
0.98mJy mas"^. RMs are shown only where the polarized flux exceeds ImJymas"-. Transverse profiles, corresponding to the labeled black bars, are shown to the 
right of each RM map (note the differences in the angular ranges). In addition to the RM profile for the view shown in each map (black), those from views with 
<I> = 0°, 45°, 90°, 135°, 180°, 225°, 270° and 315° are shown in blue for comparison (with the enclosed envelope shown in grey). In all maps and transverse 
sections the beam size is shown in the lower-left and lower-right, respectively. 
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Figure 5. RM maps at 15 GHz for jets viewed at inclinations of 10° (top) 
and 30° (bottom; cf.the upper-left panels of Figure 4; see Table 1 for a com- 
plete list of relevant parameters). For reference total flux contours are over- 
layed, in logarithmic factors of 2, with the minimum contour con'esponding 
to 0.98mJy mas"^. RMs are shown only where the polarized flux exceeds 
Imjy mas"^. To the right of each RM map, illustrative transverse RM-profiles, 
corresponding to the labeled bars are shown (black) along with those from a 
variety of azimuthal viewing angles (blue). In all maps and transverse sections 
the beam size is shown in the lower-left and lower-right, respectively. 

a factor of two in each direction (see Table 1). With this con- 
straint, the RMs are typically smaller at smaller inclinations, 
typically resulting in less depolarization. Indeed, at = 10° 
the radio core is not depolarized, though finite beam effects do 
washout the RM gradients at this location. As with our canoni- 
cal model, beyond a single beam width the transverse RM pro- 
files become indicative of the ordered, ideal resolution values. 

At larger inclinations the typical RMs increase. As in the 
canonical case, the transverse RM profiles across the polar- 



ized flux core exhibit large variations among azimuthal view- 
ing angles, and are thus unreliable measures of the circum-jet 
material. Less obvious is the fact that the intermediate trans- 
verse profile. Figure 5d, also fails to reproduce the ideal resolu- 
tion RM gradients, despite the strong correlation between RM 
profiles obtained from differing azimuthal viewing directions. 
Based upon the simple model for beam-broadened RM gradi- 
ents in Equation (25), we would anticipate a gradient roughly 
four times larger Thus the substantially reduced values of the 
RM gradient cannot be accounted for simply due to the beam- 
broadened underlying RM gradients. Nevertheless, at suffi- 
ciently large distances, where the jet is resolved in the trans- 
verse direction, the RM gradients are reproduced with high fi- 
delity, and the large-scale structure is clearly evident. 

3.1.2. Black Hole Mass 

Variations in the black hole mass serve both to assess the im- 
portance of the intrinsic astrophysical variation in the masses of 
AGNs as well as the effects of our extrapolation scheme (Sec- 
tion 2.1.1). Here we choose to keep the physical size of the re- 
gion of interest fixed, and therefore change Rmax/Raim oc M"', 
i.e., change the degree of extrapolation. Since in all cases we 
are looking at the same underlying jet structure, we are neces- 
sarily gauging the relative importance of the black hole mass, 
or more specifically the length scale it provides. 

Again we must address the problem of constructing observa- 
tionally comparable systems. This is complicated by the fact 
that m = M/Mq = 3 x 10" systems, the mass required to avoid 
any extrapolation given the physical size (15pc) and distance 
(300 Mpc) under consideration, are physically unrealistic. In 
principle, as in Figure 5 we would like to set M such that the 
total fluxes are comparable. However, this suggests the natural 
choice of fixed M, orM/Medd ocM"', which is what we choose. 
If the emission were optically thin everywhere this would give 
the desired result. In practice the total flux varies by a factor 
of 3 over variations in mass over 3 orders of magnitude, os- 
tensibly due to the subtlety different magnetic field geometry 
(see below). We might also expect the RMs to be compara- 
ble: with n cx M/r^ and B cx n\l^lj{M/r) oc M^^^/r we have 
RM - nBr oc M^l-/? and thus at fixed r, RM - M^/-. This 
is not, however, the case, with the typical RMs of less massive 
black holes being considerably larger than those of their larger 
analogs. Far from the radio core, where the beam convolved 
RMs are indicative of the true values, we find that in practice 
RM cx M"' (though, for a given mass the RMs do have the ex- 
pected dependence upon accretion rate). 

Near the black hole the toroidal and poloidal components of 
the magnetic field have similar magnitudes. Thus the results of 
Section 2.1.1 imply that at fixed physical distance /b'^ cx M. 
Since outside the jet b'' undergoes stochastic variations, we 
may therefore expect to find less ordered RM profiles at larger 
masses, exhibiting more vertical structure. As seen in Figure 
6, this is indeed the case. For m = M/Mq > 3 x 10'" vertical 
gradients can begin to dominate the RM map near the core, ex- 
tending further and becoming more prominent as m increases. 
This vertical structure is an intrinsic effect, appearing in the 
ideal resolution RM maps as well. This is in contrast to RM- 
gradient reversals near the core for m <3 x 10'", which are 
exclusively due to finite beam effects. 

The transverse gradients near the radio core at lower masses 
are clearly not indicative of the true, nearly linear RM gradi- 
ents. However, for ;7i > 3 x 10'" the transverse RM profiles in 
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Figure 6. RM maps at 15GHz for jets produced by black holes with a variety of masses {m = M/Mq is listed in the lower-right of each RM map; see Table 1 for a 
complete Ust of relevant parameters). Of particular note is the fact that at 3 X 1O"M0 no extrapolation is required. Our canonical model is shown in the lower-left. 
For reference total flux contours are overlayed, in logarithmic factors of 2, with the minimum contour corresponding to 0.98 mjy mas"-. RMs are shown only where 
the polarized flux exceeds ImJymas"^. To the right of each RM map, illustrative transverse RM-profiles, corresponding to the labeled bars are shown (black) along 
with those from a variety of azimuthal viewing angles (blue). In all maps and transverse sections the beam size is shown in the lower-left and lower-right, respectively. 



14 



BRODERICK & MCKINNEY 



Figure 6c & 6f are more robust, exhibiting a clearly identifi- 
able gradients. While the RM gradients along the nearby sec- 
tions are well below the true values, they are consistent with the 
beam-broadened values anticipated by Equation (25). Since in 
practice the jet width is not known a priori, determining the 
true gradients along these sections is not possible for AGN jets 
associated with the most massive black holes expected in na- 
ture. Nevertheless, in this case the presence of a robust gradi- 
ent, even if unresolved, is indicative of ordered toroidal fields. 

3.1.3. Magnetic Dissipation in the Faraday Screen 

Finally, we consider the dependence of the RMs upon the 
rate at which magnetic dissipation heats the circum-jet mate- 
rial. The transition from hot ikT > m^c-) gas within the jet to 
cool (kT < nieC') gas surrounding the jet is typically rapid, oc- 
curring over less than a jet width. Nevertheless, since f{T) is 
such a strong function of kT /m^c- (see Equation (21)), the lo- 
cation of this transition can substantially affect the typical RMs 
produced near the jet. Since b~ /%Tipc- > 1 within the jet, effi- 
cient dissipation of magnetic energy in the jet sheath can easily 
heat the circum-jet material sufficient to suppress Faraday ro- 
tation. 

This is explicitly seen in Figure 7, which shows the efficient 
heating model for a number of maximum heating rates e (re- 
call, we set the maximum internal energy to eb~/d,TT). When 
£=1.0 the typical RMs are reduced by more than two orders of 
magnitude. In contrast, setting e < 0.01 results in nearly iden- 
tical values to those from our canonical model. In all cases, 
with the exception of regions near the core, the basic structure 
of the RM maps and their transverse profiles are qualitatively 
similar. 

As for the canonical case (shown in the lower-right panel of 
Figure 7), where the transverse structure of the jet is resolved 
the RMs are indicative of the true values. Similarly, the trans- 
verse profiles near the core are typically unreliable. An excep- 
tion is the appearance of a robust gradient at this location when 
e = 0.1 (Figure 7f). In this case the resulting beam-convolved 
gradients are indeed consistent with those expected from Equa- 
tion (25). In contrast, despite exhibiting a general trend, the 
RMs in Figure 7i are roughly an order of magnitude lower than 
would be expected. As a result, it may be difficult to assess 
the fidelity of an observed transverse RM gradient when the jet 
width is insufficiently unresolved. 

3.2. Dependence of the RM upon the Observation Parameters 

In addition to the parameters of the AGN model, the mea- 
sured RM can depend upon how they are observed as well. 
Here we address two ways in which this occurs: observing fre- 
quency and beam size. In order to isolate the effects of each, we 
treat them separately by fixing one while varying the other All 
of these are essentially finite beam effects; in all cases shown 
in Figures 9 & 8 the intrinsic, ideal resolution RM maps are 
identical. 

3.2.1. Beam Size 

We have alluded to the importance of beam size a number 
of times already; generally finding that where the beam does 
not resolve the transverse structure of the jet the transverse RM 
profiles are inaccurate. Here and in Figure 8 we make these 
statements explicit by convolving the same image (correspond- 
ing to our canonical model; see Table 1 for details) with a num- 
ber of different beams. These correspond to the highest typi- 
cal VLBA resolution at 5 GHz (0.9 mas), highest typical VLBA 



resolution at 43 GHz (0.1 mas), highest design resolution for 
VLBI using the upcoming space-based radio telescope VSOP- 
2 (43 GHz, 0.04 mas, Tsuboi et al. 2009) and the computational 
resolution (~ 0.02 mas). These can also be compared with our 
canonical model (e.g. Figure 4), which has a beam typical of 
VLBA observations at 15 GHz. 

The 0.9 mas and 0.3 mas beams are typical of matched- 
resolution RM studies, and thus provide insight into the ability 
of the existing literature to constrain the parameters of simu- 
lated jets. VSOP-2 promises to achieve beam sizes of 0. 1 mas 
at 8 GHz, roughly a factor of two improvement over the VLBA. 
Thus, 0.1 mas beam is indicative of what may be possible in 
the near future. Finally, VSOP-2's maximum resolution gives 
a lower-limit on realistically attainable beam widths over the 
next decade, though in this case the observed RMs would be 
obtained only over a relatively small wavelength range. 

As previously stated, the ideal resolution case produces clear, 
nearly linear transverse gradients along the length of the jet. 
VSOP-2 should be able to probe gradients as close as 0.3 mas 
to the radio core. Similarly, high frequency observations at the 
VLBA faithfully reproduce the transverse RM gradients in the 
optically thin region, though near the radio core large spurious 
fluctuations develop. 

In contrast, for large beams, while apparently linear gradi- 
ents exist, these are in some cases clearly anomalous, arising 
from the breadth of the beam instead of the structure of the 
jet. This is most clearly evident in the bottom two transverse 
profiles. Figure 8k & 81, which show profiles with a number 
of different gradients, in some cases with the wrong sign. In 
Figure 8k, the presence of inverted gradients provides a strong 
case against interpreting these as due to the true structure in the 
RM map. However, the average profile in this case is roughly 
consistent (within 50%) with that anticipated by Equation (25). 
In contrast. Figure 81 exhibits RM gradients which are all of 
the correct sign, though also all well below the values expected 
even after accounting for the beam-broadening of the profiles. 
Far from the radio core, despite being unresolved, the trans- 
verse gradients are strongly correlated in azimuthal viewing 
angle and are well-described Equation (25), producing gradi- 
ents roughly a factor of 4 smaller than the true values. Thus, 
despite the inability to infer the actual RM values due to the 
failure to resolve the jet transverse structure, the structure in 
the RM map is indeed indicative of the ideal-resolution RMs 
more than a beam width from the radio core. 

3.2.2. Frequency 

While we have defined the RM in the differential sense, 
where the A^-law holds in the ideal-resolution RM maps, we 
might anticipate that frequency effects may be neglected. How- 
ever, as seen in Figure 9 (cf. 4), near the radio core this is not 
the case generally. Despite having identical intrinsic RMs, dif- 
ferences in the photosphere size and beam depolarization as- 
sociated with different observing frequencies can dramatically 
alter the beam-convolved RM maps. Note that here we fix the 
beam size to 0.3 mas despite the change in observing frequency. 

Firstly, the increase in beam depolarization at longer wave- 
lengths severely restricts the region over which the RMs may 
be measured. At 5 GHz, the jet is a few beams across, and 
begins to suffer from inadequate transverse resolution. As a 
consequence, even where the jet is optically thin, the beam- 
convolved RM map exhibits substantial variations about the 
true values. The RM profile shown in Figure 9b is roughly 
a factor of two below the true values, though in agreement with 
those predicted by Equation (25). By 8 GHz the jet is suffi- 
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Figure 7. RM maps at 15 GHz for jets in wliich the jet slieatli is efficiently heated via magnetic dissipation (e is listed in the lower-right of each RM map; see 
Table 1 for a complete list of relevant parameters). For each the internal energy is extrapolated using a„ = -1.3 until it reaches efo^/S-Tr (generally, this occurs quite 
close to the black hole). The upper-left, upper-right and lower-left panels correspond to £ = 1.0, 0.1 and 0.01, respectively. For comparison the canonical model 
{flu = —8/3) is shown in the lower-right panel. For reference total flux contours are overlayed, in logarithmic factors of 2, with the minimum contour corresponding to 
0.98mJy mas"^. RMs are shown only where the polarized flux exceeds ImJymas"-. To the right of each RM map, illustrative transverse RM-profiles, con'esponding 
to the labeled bars are shown (black) along with those from a variety of azimuthal viewing angles (blue). In all maps and transverse sections the beam size is shown 
in the lower-left and lower-right, respectively. 
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Figure 8. RM maps of our canonical model at 15 GHz as observed with a number of different beam sizes (beam FWHMs are listed in the lower-right of each 
RM map; cf. the upper-left panel of Figure 4; see Table 1 for a complete list of relevant parameters). In particular, we show the map at the computed pixel 
resolution (0.02 mas, upper-left), the highest design resolution of VSOP-2 (0.04 mas, upper-right), the highest typical VLBA resolution (0.1 mas, lower-left) and the 
typical VLBA resolution at 5 GHz (0.9 mas, lower-right). For reference total flux contours are overlayed, in logarithmic factors of 2, with the minimum contour 
corresponding to 0.98mJy mas"-. RMs are shown only where the polarized flux exceeds ImJymas"^. To the right of each RM map, illustrative transverse RM- 
profiles, coiTesponding to the labeled bars are shown (black) along with those from a variety of azimuthal viewing angles (blue). In all maps and transverse sections 
the beam size is shown in the lower-left and lower-right, respectively. 
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Figure 9. RM maps of our canonical model as observed v = 5 GHz, 8 GHz, 22 GHz and 43 GHz (observation frequecies are listed in the lower-right of each RM 
map; cf. the upper-left panel of Figure 4; see Table 1 for particular parameters). In all cases the beam size is fixed to 0.3 mas. The beam-convolved spectral index 
map is shown in the upper-left inset of each RM map. in which the colors range from —0.5 (dark blue) to 0.5 (dark red), indicating the size of the photosphere at each 
frequency. For reference total flux contours are overlayed. in logarithmic factors of 2 (factors of 8 in the a map insets), with the minimum contour con'esponding to 
0.98mJy mas"^. RMs are shown only where the polarized flux exceeds ImJymas"^. To the right of each RM map. illustrative transverse RM-profiles, con'esponding 
to the labeled bars are shown (black) along with those from a variety of azimuthal viewing angles (blue). In all maps and transverse sections the beam size is shown 
in the lower-left and lower-right, respectively. 
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ciently resolved (given our polarized flux limit of 1 mJymas"-) 
that the gradients in the optically thin portions of the jet (e.g. 
the top two transverse profiles) reproduce the true gradients 
with high fidelity. 

Secondly, the size of the radio photosphere shrinks as fre- 
quency increases, with a corresponding decrease in the conse- 
quences of the rapid variations in the spectral index. At 5 GHz 
the photosphere extends roughly 0.5mas, larger than the beam 
itself, while at 43 GHz it is comparable to the resolution with 
which we computed the intrinsic RM maps (0.02 mas). The de- 
crease in the importance of the optically thick core is seen most 
dramatically in the appearance of the consistent transverse gra- 
dients close to the core. Figure 91, which are roughly an order 
of magnitude below the true values, consistent with Equation 
(25). 

When the beam size is increased to 0.9 mas all of the profiles 
at 5 GHz fail to reproduce even the beam-broadened RM gra- 
dients. Nevertheless, above 8 GHz, as in Figure 8j, transverse 
profiles far from the core exhibit robust gradients indicative of 
the expected values for transversely unresolved jets. This is 
also true at intermediate distances from the radio core (Figure 
9h & 9k) for frequencies above 22 GHz. Finally, at 43 GHz, the 
highest frequency we consider, even the core RM profiles are 
well-defined and in agreement with the values anticipated by 
Equation (25), though roughly two order of magnitudes smaller 
than the true values. Thus, for the purposes of detecting robust 
RM gradients associated with poorly resolved structure within 
the Faraday screen, pushing to high frequencies alone can par- 
tially ameliorate resolution issues (though, as addressed explic- 
itly in Section 3.4, it may be difficult to distinguish between 
nearby and distant Faraday screens in this case). 

The differences between the RM maps near the radio core 
due to beam convolution necessarily implies that the A^-law 
is violated in this region. On the other hand, when the trans- 
verse jet structure is sufficiently resolved and the observed 
RMs agree with the intrinsic RMs, as occurs far from the core 
in the optically thin region, then similarities between the RMs 
implies that the A--law holds. These regions also differ in how 
well the beam-convolved RM maps reproduce the true values. 
Hence, this suggests that at least for our simulated AGN jets, 
the A^-law provides a rough measure of the accuracy of the RM 
gradients when the RMs are sufficiently resolved^. 

3.3. Dependence of the RM upon the Jet Structure 

The primary motivation for studying the RM distributions of 
AGN jets has been to infer the structure and parameters of the 
Faraday screen. In the context of the 3D simulation, neither of 
these are free parameters, and we have obtained a characteris- 
tic, bilateral structure for the RM maps, presumably as a conse- 
quence of large-scale structure within the Faraday screen. Here 
we address how these morphological features are related to the 
structural properties of the Faraday screen itself. We do this 
by artificially removing components of the magnetic field & 
velocity and recomputing the RM map, while using our canon- 
ical model to produce the jet emission. 

3.3.1. Magnetic Field Structure 

The relationship between the poloidal and toroidal magnetic 
field has already been indirectly demonstrated to have a signif- 

' This conclusion in specific to our jet models and may not be generic. For 
example, unresolved point sources are capable of producing RMs which appear 
to follow the A- though are erroneous. This does not appear, however, to be 
the case for the highly-ordered structures that appear within our simulated jets. 
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Figure 10. RM maps at 15 GHz when toroidal component of the comoving 
magnetic field (top) and poloidal component of the velocity (bottom) of the 
plasma within the Faraday screen is made to artificially vanish (cf. the upper- 
left panel of Figure 4; see Table 1 for a complete list of relevant parameters). 
For reference total flux contours are overlayed, in logarithmic factors of 2, with 
the minimum contour corresponding to 0.98mJymas"^. RMs are shown only 
where the polarized flux exceeds ImJymas"^. To the right of each RM map, 
illustrative transverse RM-profiles, corresponding to the labeled bars are shown 
(black) along with those from a variety of azimuthal viewing angles (blue). In 
all maps and transverse sections the beam size is shown in the lower-left and 
lower-right, respectively. 

icant effect upon the structure of the RM map. Large poloidal 
fields (e.g., those associated with m = 3 x 10") can result in the 
loss of the transverse symmetry apparent outside of the radio 
core. This is illustrated explicitly in the top panel of Figure 10, 
which presents the RM map when the toroidal component of 
the magnetic field in our canonical model is artificially made 
to vanish. Significant transverse gradients vanish completely 
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everywhere within the jet (where it is sufficiently resolved), 
confirming that nearly-linear, resolved transverse gradients are 
due to large-scale, ordered toroidal magnetic fields within the 
Faraday screen. 

The remnant RM is that associated with the poloidal field, 
and is thus roughly 3 orders of magnitude smaller than that 
in our canonical model (see. Figure 4). This is a function of 
the black hole mass, with jh^ cx M"' at fixed physical dis- 
tance, and thus the perturbations due to poloidal fields will be 
larger for more massive systems. Nevertheless, it appears that 
for typical blazars and BL Lac objects the fixed offset in the 
transverse profile due to the presence of poloidal fields will be 
dominated by the uncertainties in the location of the jet spine 
(not necessarily located along the jet axis), viewing angle and 
finite beam effects for single RM profiles. It remains to be 
seen if it is possible to reduce these uncertainties by combining 
multiple profiles, making use of the strong correlation in these 
quantities along the jet axis. 

33.2. Velocity Structure 

The large Lorentz factors within AGN jets suggest that the 
velocity structure of the Faraday screen may also be important. 
The degree to which this occurs in practice depends upon the 
location of the Faraday screen and how coupled it is to the jet 
itself. 

Setting the poloidal velocity to vanish, shown in the bottom 
panel of Figure 10, made little difference to the overall mor- 
phology of the RM map. It does, however, decrease the typical 
RMs by a factor of roughly 4. This is not unexpected if the 
typical Lorentz factors within the Faraday screen are on the or- 
der of 2 (a fact that will be borne out in the following section); 
from Equation (21), 
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which when O < Vl/T reduces to r~~. Unlike the canonical 
case, when the poloidal velocity is made to vanish the trans- 
verse RM profiles exhibit a well-defined, nearly linear gradi- 
ents near the core. These are roughly an order of magnitude 
smaller than the true RMs, though in agreement with those ex- 
pected by Equation (25). 

In contrast, the toroidal velocities appear to have little effect 
beyond the radio core. This is not particularly surprising since 

from Section 2. LI we have that (3'^ /P'' ^ (GM/c~r)~ , and 
thus at the large distances (in terms of the gravitational radius 
of the black hole) seen here we anticipate the toroidal velocities 
to be roughly 5 orders of magnitude smaller than the poloidal 
values within the screen. Thus the kinds of effects associated 
with large-scale relativistic motions described in Broderick & 
Loeb (2009) are likely to be unimportant in the case of isolated, 
roughly stable jets. 

3.4. Location and Structure of the Jet Faraday Screen 

While we have described how some of the most commonly 
discussed features of the circum-jet environment imprint them- 
selves upon the RM map, we have said little about the location 
and structure of the Faraday rotating medium. Here we address 
this, identifying what portion of the outflow is responsible for 
generating the Faraday rotation and how it is related to the jet 
itself. In addition, we address how this may be distinguished 
from non-local models of the Faraday screen. 
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Figure 11. Plasma quantities along an example line of sight, in comparison to 
the location of the Faraday screen. The locations of the jet, defined by pc? < 
ub = fc^/S-TT, and its surrounding sheath, defined by pg < ub/3, are represented 
by the dark and light gray shaded regions, respectively. This corresponds a 
Une of sight responsible for producing the RM at an angular' position of AQ^ = 
Imas and Af2y = 1.5 mas, and is the second from the bottom Une of sight 
shown in Figure 2, with the intensity (soUd black line) and rotation measure 
(dashed, blue Une) shown explicitly in the bottom panel. 

A number of relevant plasma quantities, in addition to and 
the RM, along a representative line of sight are shown in Fig- 
ure 1 1 , including the poloidal and toroidal components of f3 
and b, the Lorentz factor and temperature. Generally, the Fara- 
day screen is quite close to the jet, lying within two jet widths 
of the jet axis. The bulk of the Faraday rotation occurs within 
the strongly magnetized jet sheath, defined here by when the 
gas pressure is less than 1 /3 of the magnetic pressure (i.e., 
the plasma /? < 1/3), though at larger inclinations a significant 
fraction may occur within the surrounding wind. While the 
Faraday screen and emission region are well-separated, they 
are both parts of a larger MHD structure. As illustrated in Fig- 
ure 11, the jet sheath is a smooth extension of the jet itself; 
i.e., there is no well-defined spine-sheath structure despite our 
identification of a jet "sheath" (cf. Laing & Bridle 2002b, a). 

The entire rotating region is unbound in the particle, thermal 
fl«(i magneto-thermal senses, i.e.. 



2GM\ 
c^r I 



r>i, 



p 



pu 



\ c r J p 
2GM\ {p + ug+pg) r^P''-l/f3-h 



(28) 



c^r J 



respectively, where T,'' is the contravariant poloidal, covariant 
time component of the ideal MHD stress-energy tensor, u, is 
the covariant time component of the 4-velocity, and u'' is the 
contravariant poloidal component of the 4-velocity. Therefore, 
the Faraday screen is clearly associated with the large-scale 
outflow. Despite this, the Lorentz factor of the Faraday rotat- 
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ing material is generally much lower than the jet interior, F ~ 2, 
due in part to the strong dependence upon F of the integrand in 
Equation (6). 

The RM is dominated by contributions close to the jet and 
not by the disk wind. In the simulations, the disk wind is not as 
well-collimated as the jet due to the finite extent of the accre- 
tion flow. While the highly relativistic jet becomes mostly bal- 
listic when losing lateral support due to the finite disk extent, 
the weakly relativistic disk wind fills the space as a more spher- 
ical outflow. Even the part of the disk wind launched from near 
the black hole has a significant loss of toroidal field strength 
with radius compared to the relativistic jet due to the overall 
expansion of the disk wind. The simulations suggest that a 
disk extent of order pc-scales may be required to have a well- 
collimated disk wind that significantly contributed to the RM 
on pc-scales. 

Throughout the duration of the 3D GRMHD simulation we 
have used, there is no evidence for the presence of reversed 
field polarity as in the "magnetic tower" picture discussed by 
Mahmud et al. (2009). They suggested that the time-dependent 
sign of the RM in BL Lac object B1803H-784 can be explained 
best by magnetic tower models containing nested, reversed he- 
lical fields, in which the inner field changes (in time) its degree 
of winding relative to the outer return field. In the GRMHD 
simulations, at no point along the line of sight does the toroidal 
component of the comoving or lab-frame magnetic field re- 
verse direction It appears unlikely that a significant fraction 
of the observed RM can be due to a distant, reversed region, 
as in a disk wind. Furthermore, such reversed magnetic tower 
geometries appear to be at odds with jets launched by accre- 
tion disks, for which the interior of the jet is universally more 
tightly wound due to the monotonically decreasing radial Ke- 
plerian rotation profile. One might be able to invoke a slowly 
rotating black hole connected to time-varying launching points 
in the Keplerian disk. However, this requires a slowly rotat- 
ing black hole with a/M < 0.4, at which the total power of the 
disk wind may dominate that of the black hole jet (McKinney 
2005) leaving the emission dominated by a quasi-spherical re- 
gion instead of a collimated jet. Sign changes could potentially 
be due to changes in the overall polarity of a dipolar field near 
the black hole, although this may be at odds with the long-term 
constancy of the sign of the RM gradient. On the other hand, 
the presence of instantaneous erroneous RM gradients (even 
sign changes) due to the finite beam size can produce qualita- 
tively similar behaviors near the radio core, and might be re- 
sponsible for the apparent time-dependence, though we should 
note that in B1803H-784 the reversal occurred far from the ra- 
dio core. An investigation of RMs that considers these issues, 
such as considering RMs from time-dependent simulations, is 
left for future work. 

The proximity of the Faraday screen to the jet makes a strong 
prediction regarding the morphology of the RM maps. The 
presence of a large-scale radially correlated, toroidally domi- 
nated magnetic field within the jet sheath implies that the RM 
map will be strongly correlated along the jet. This has dra- 
matic consequences for the sgn(RM) maps, resulting in a very 
particular bilateral structure. 

Example sgn(RM) maps are presented for a variety of jet 
model and observational parameters in Figure 12. For two 
different beam widths, chosen to represent resolved and un- 
resolved jets, (0.3 mas and 0.9 mas on the top and bottom, re- 
spectively) we present the sgn(RM) maps representative of the 
various parameters we have considered in Sections 3.1-3.3. 
Specifically, for each beam size we show the sgn(RM) asso- 



ciated with our canonical model (upper left), at 43 GHz (lower 
left), high and low inclinations (upper and lower center-left, 
respectively), a black hole with m = 3 x 10'" (upper center), 
an efficiently heated Faraday screen (lower center), vanish- 
ing toroidal magnetic field (upper center-right) and vanishing 
poloidal velocity (lower center-right). The sgn(RM) morphol- 
ogy is insensitive to the properties of the poloidal magnetic 
field and the toroidal velocity, though both can alter the details 
of the sgn(RM) map. 

With the exception of a vanishing toroidal magnetic field, 
all of the jet models listed above exhibit the afore mentioned 
bilateral morphology throughout the optically thin portion of 
the jet (as with the RM profiles, near the radio core departures 
from this structure can occur for some azimuthal orientations). 
This suggests that if the Faraday screen is close to the jet, the 
observation of such a structure in the RMs of the optically thin 
portions of AGN jets is a robust indicator of the presence of 
ordered, toroidal magnetic fields. However, the strength of this 
evidence does depend upon the particulars of the jet model, and 
more importantly the beam resolution. 

At small inclinations, insufficient radial range may exist to 
argue unambiguously for the correlated bilateral morphology. 
At large inclinations, the transverse structure of the jet may re- 
main unresolved, leading to the previously-mentioned issues 
surrounding beam convolution. This is illustrated in the large- 
beam sgn(RM) maps (Figure 12, bottom), in which the unre- 
solved maps show relatively little structure in all cases. While it 
still appears possible to determine if ordered toroidal magnetic 
fields exist within the Faraday screen, the large beam limits the 
amount of asymmetric structure present when toroidal fields 
are sub-dominant. This difficulty vanishes when the transverse 
jet structure is resolved. 

The morphology of the sgn(RM) maps may also be used to 
observationally probe the location of the Faraday screen, ruling 
out contributions from distant, randomly oriented foreground 
clouds. While pc-scale clouds have been detected via free-free 
absorption in a handful of sources (Jones et al. 1996; Walker 
et al. 2000), typically they are incapable of reproducing the 
large-scale radial correlations associated with Faraday rotation 
within the jet sheath. This may be argued on general grounds: 
Any model for the foreground Faraday screen must have struc- 
ture on scales comparable to the jet width in order to mimic the 
observed transverse gradients. However, since in this case the 
Faraday screen is not associated with the jet itself, this implies 
that there must be structure parallel to the jet as well. Thus, we 
do not expect to see the bilateral morphology associated with 
the Faraday rotation due to a jet sheath. 

To illustrate this explicitly, in the right-most panels of Fig- 
ure 12 we show two examples of random distant foreground 
screens, in which dRM./dl is given by a Gaussian random 
field with power-law power spectra, Pk oc k"' . To emphasize 
the diagnostic ability of the sgn(RM) map, we subtract off the 
flux-weighted average RM, which is degenerate with a Galac- 
tic contribution, in order to ensure that the sgn(RM) map ex- 
hibits significant structure. These realizations are not intended 
to represent physical models for the plasma conditions within 
the Faraday screen. Rather, they simply provide two well- 
defined examples of the beam-convolved sgn(RM) maps asso- 
ciated with foreground screens that are uncorrected with the jet 
itself. In keeping with this motivation, we choose a flat spec- 
trum (fl, = 0), which necessarily produces considerable struc- 
ture on angular scales comparable to the beam, and a red spec- 
trum {ur = -3) which typically generates structure on the largest 
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Figure 12. sgn(RM) for an illustrative subset of the cases discussed in the text, compared to two examples of random foreground screens (rightmost maps). These 
are shown for two beam sizes, con'esponding to resolved (top) and unresolved (bottom) jets. See the text and Table 1 for a complete list of relevant parameters for 
each, though the relevant defining parameter is listed as it appears in Figures 5-10 in the lower-right corner of each sgn(RM) map. For reference total flux contours 
are overlayed. in logarithmic factors of 2, with the minimum contour corresponding to 0.98mjymas"^. RMs are shown only where the polarized flux exceeds 
ImJy mas"^. In all maps the beam size is shown in the lower-left comer. 
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relevant scale (e.g., the jet length). 

The flat-spectrum foreground screen produces sgn(RM) 
maps that are similar in structure to those arising when the 
toroidal magnetic field within the jet sheath vanishes. As a con- 
sequence, our discussion of the sgn(RM) maps as indicators of 
the presence of ordered toroidal fields applies to the distinc- 
tion between local and distant, flat-spectrum Faraday screens 
as well. That is, the presence of ordered, bilateral structures in 
the optically thin portions of the jet necessarily excludes this 
kind of distant screen. 

The red-spectrum foreground screen is more problematic. In 
this case large-scale structures characteristic of red spectra are 
apparent, frequently extending substantial distances along the 
jet axis. As a result, distinguishing between red-spectra fore- 
grounds and nearby Faraday rotation is difficult at large beam 
widths. This situation is ameliorated somewhat by smaller 
beams, at which not only the general morphology but also the 
detailed structure of the sgn(RM) maps are apparent (cf. the 
canonical and red-spectrum models at small beams in Figure 
12). Despite this, if the jet sheath is efficiently heated, this com- 
plication may remain, making the distinction between fore- 
ground and nearby Faraday rotation for any individual source 
challenging. 

3.5. Implications for Jet Models 

Because the RM depends upon the thermal particles that are 
directly modeled in jet theory, it provides a more robust han- 
dle upon the parameters of jet models. Thus, in addition to 
addressing prior suggestions about the source of the RMs asso- 
ciated with AGN jets, we may also confront the observational 
applicability of our present understanding of how relativistic 
jets are produced and propagate. 

First among these is the relationship between the nonthermal 
particles and the other jet parameters. By normalizing the ac- 
cretion rate via the typical RM observed (~ 10^'radm"-), we 
have been able to determine the efficiency with which the elec- 
tromagnetic energy within the jet is dissipated into nonthermal 
particles producing the emission. For our canonical model we 
find that rj = 0.02 reproduces typical observed fluxes, meaning 
that the nonthermal particles constitute approximately 2% of 
the energy density within the jet. This is roughly an order of 
magnitude larger than internal energy of the thermal particles 
within the jet. Therefore, when combined with RM measure- 
ments, their luminosity implies that at pc-scales AGN jets are 
Poynting-dominated outflows of nonthermal, relativistic parti- 
cles. This is consistent with some studies of kpc-scale AGN 
jets, finding that even at those scales large-scale fields persist 
and dominate the plasma thermal energy (see, e.g., Lapenta & 
Kronberg 2005b). 

Second, constraining the plasma parameters of the Faraday 
screen via RM profiles is complicated and will generally re- 
quire careful modeling of the screen properties, e.g., via simu- 
lations like that we have presented here. This is because there 
are significant degeneracies between a number of otherwise 
poorly constrained properties; namely the jet orientation, bulk 
velocity of the rotating medium, thermal plasma density and 
temperature and the strength of the magnetic field. As seen 
in Figure 11, many of these are varying substantially through- 
out the screen. Furthermore, simplifying assumptions, such as 
equipartition, can be wrong by orders of magnitude. While it 
may be tempting to calibrate simple, qualitative models using 
simulations like the one we present, the relationship between 
the circum-jet material and the fast jet core is a strong func- 
tion of black hole spin and other model parameters. Thus, 



efforts to quantitatively infer the properties of the circum-jet 
material will require additional, self-consistent jet formation 
simulations. 

Third, we nevertheless find that the ordered toroidal mag- 
netic field in the jet and its sheath were crucial to producing 
emission and Faraday rotation distributions that comport with 
observations. This is consistent with simplified, analytical ef- 
forts to explain these properties with a strongly-magnetized jet 
containing large-scale ordered toroidally dominated magnetic 
fields (see, e.g. Lyutikov et al. 2005). If AGN jets correspond 
well to GRMHD jet models like that presented here, observa- 
tions of the RM gradients coupled with constraints upon the jet 
orientation, particle acceleration and heating rates within the jet 
and M as estimated by the disk thermal spectra and luminosity 
provide a means to estimate the jet power per unit mass accre- 
tion rate. For jets powered via the Blandford-Znajek mecha- 
nism, this quantity is a strong function of the black hole spin 
(McKinney 2005; Tchekhovskoy et al. 2010), and hence may 
provide a means to measure the black hole spin magnitude for 
a large number of AGNs. 

Fourth, it has been suggested that variations in the polariza- 
tion angle across the jet indicate changes in the direction of the 
intrinsic magnetic field. Specifically, in a number of sources 
(e.g., 3C 272 and 1055H-018) the polarization is aligned with 
the jet near the axis and becomes orthogonal to the jet at the 
edges (Attridge et al. 1999; Pushkarev et al. 2005; Zavala & 
Taylor 2005; Gomez et al. 2008). Often this is considered 
evidence for a toroidally dominated jet core and a poloidally 
dominated jet sheath. However, studies of purely toroidal jets 
have qualitatively reproduced these polarization distributions 
as a consequence of projection effects and relativistic aberra- 
tion (Lyutikov et al. 2005; Zakamska et al. 2008). Seen ex- 
plicitly in Figure 3, our work also qualitatively reproduces the 
observed polarization structure without the need to appeal to 
variations in the intrinsic magnetic field geometry. 

Finally, while the jet simulation is fully 3D, and this non- 
axisymmetric structure is clearly evident in the jet flux dis- 
tributions, the resulting RM profiles are remarkably uniform 
(where they are sufficiently resolved). Thus, despite substan- 
tial departures from axisymmetry observed in the jet core, the 
assumption of axisymmetry in the surrounding medium seems 
well-justified. 

4. CONCLUSIONS 

It has now become possible to quantitatively compare pc- 
scale observations of AGN jets with self-consistent fully 3D 
numerical GRMHD simulations of jets. Due to the need to 
model the nonthermal particles, direct comparisons between 
the emission and polarization properties of simulated and real 
jets are not possible, generally requiring additional post hoc as- 
sumptions about the poorly understood process of particle ac- 
celeration within the jet. In contrast, the RMs associated with 
the circum-jet outflows depend only upon the simulated quan- 
tities. For these reasons AGN jet RMs provide a potentially 
more robust diagnostic of the simulation parameters. 

Even with relatively low accretion rates, ^ 10~^MEdd, we are 
able to reproduce the typical RMs of lO^radm"^ observed in 
pc-scale jets. These are generally proportional to Ivp'l", and 
thus even moderate increases in the accretion rate results in 
substantial increases in the typical RMs. As a consequence 
it is possible to produce the majority of the observed Faraday 
rotation via propagation through the circum-jet material alone. 
Were the bulk of the observed RMs to be produced far from the 
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jet, this would require either Ms well below that necessary to 
produce the observed emission or very efficient dissipation of 
the magnetic field immediately outside of the jet (and thus very 
hot circum-jet plasma). 

While well-separated from the emission region, the bulk of 
the Faraday rotation occurs in the immediate vicinity of the jet, 
within the surrounding jet sheath. This consists of a highly 
magnetized, relativistic outflow, which though moving more 
slowly than the jet core is more properly thought of as a con- 
tinuous extension of the jet itself than a separate portion of the 
outflow. As a consequence, like the jet core, the magnetic field 
in the Faraday rotating jet sheath is dominated by a large-scale, 
ordered toroidal magnetic field. 

The RM maps associated with our canonical model ex- 
hibit well-defined transverse RM gradients, associated primar- 
ily with the ordered toroidal magnetic field. Thus, our re- 
sults suggest that efforts to probe the properties of toroidal 
magnetic fields in the circum-jet regions are consistent with 
existing 3D jet formation simulations, which self-consistently 
model the fast jet core and the surrounding Faraday rotating 
material. Furthermore, while some variation in the precise gra- 
dients across the jet exist due to non-axisymmetric jet struc- 
tures, these are generally small in comparison to the gradients 
themselves, implying that RM observations can produce well- 
defined estimations of the toroidal field strength and structure. 

The RM map details are also sensitive to the poloidal ve- 
locities, though in a way that is roughly degenerate with the 
dependence upon the mass accretion rate. However, at pc dis- 
tances the comoving poloidal magnetic fields and toroidal ve- 
locities are sufficiently small that neither produce perturbations 
in the transverse RM profiles that exceed the variations associ- 
ated with the non-axisymmetric structure of the jet. While this 
is dependent upon black hole mass, for reasonable values this 
implies that efforts to probe the poloidal magnetic field and 
toroidal velocities may be challenging. 

Accurately measuring the RMs of AGN jets depends criti- 
cally upon the ability to resolve the jet's transverse structure. 
This is generally violated in the radio core, and thus core RMs 
and RM-gradients may be typically unreliable. When the beam 
is moderately larger than the jet width finite beam effects re- 
sult in a suppression of the RM relative to the true values, 
though well-defined gradients may still exist. When the beam 
is much larger, or when there is significant spectral or polar- 
ization structure on sub-beam scales, the resulting RMs can 
be dramatically different from the true values. This motivates 
pushing towards higher frequencies and higher resolutions, as 
well as focusing upon resolved regions more than a beamwidth 
away from the radio core. For these purposes, space-based 
VLBI experiments like VSOP-2 are likely to be particularly 
useful. 

The large-scale correlations along the jet indicative of or- 
dered toroidal fields are present even in the sgn(RM) maps, 
producing a distinctive, bilateral morphology. Models which 
invoke Faraday rotation far from the jet generally have diffi- 
culty reproducing this structure, and may therefore be excluded 
by sgn(RM) measurements. Despite the fact that sufficiently 
red fluctuations in the foreground RMs can mimic the long- 
range order indicative of toroidal fields, sufficiently resolved 
RM observations may robustly verify the presence of Faraday 
rotation in the vicinity of the jet. 

While we have made an attempt to explore the various pa- 
rameters critical to computing the RMs of simulated jets, by 
no means has our investigation been exhaustive. The most sig- 
nificant limitation is the restriction to a single time step of a sin- 



gle jet simulation with a single prescription for the nonthermal 
particles. As a consequence we are not able to address the un- 
certainties in the RMs associated with the substantial variabil- 
ity in the accretion onto the black hole. More importantly, we 
cannot address how interactions with the surrounding medium, 
known to be important in many cases, will affect our results. 
Where these interactions significantly alter the properties of the 
jet they should also make substantial changes to the particulars 
of the RM maps. However, since the general structure of the 
RMs arise from a direct relationship with the order inherent 
in the jet's magnetic field, we would expect the generic mor- 
phology of the RMs to persist as long as the jet is toroidally 
dominated. 
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APPENDIX 
RELATIVISTIC ABERRATION 

In Section 2.3.2 we described how we performed the po- 
larized radiative transfer through relativistically moving bulk 
flows. Here we present how the polarization direction changes 
due to relativistic aberration, and thus specifying the orienta- 
tion of the parallely propagated orthonormal tetrad used to de- 
fine the Stokes parameters. We will make use of four-vectors 
and covariant analysis, though the result will be quoted in terms 
of the three vectors we have used to define the relevant plasma 
quantities. 

We choose to orient the Stokes parameters so that the pro- 
jected, positive z-axis is "up" in the image plane. This is neces- 
sarily orthogonal to the wave four- vector, k" = (l, k), and thus, 

e" = (0,e)-Cr, e = z-(z-k)k, (Al) 

in which C is an arbitrary scalar Generally, e" is only defined 
up to C, corresponding to the ambiguity associated with the 
time-component of the polarization four-vector in any given 
frame*^. Thus we are free to require e" be purely space-like in 
the jet frame, i.e., with u" = (F, Tf3), we set C such that u"ea = 0, 
obtaining, 

(A2) 

l-/3-k 

This now defines the "up" direction, as seen in the jet frame. 

In the emitter's frame, the natural orientation for the Stokes 
parameters is aligned with the projected magnetic field, which 

is given by b -(b • k )k , where all primed quantities are in the 
jet frame (though note that k' is the spatial components of k" in 

^ Alternatively, this may be thought of as the gauge freedom associated with 
the four-vector potential. 
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the jet frame, and thus is not unit-normaUzed). Thus, 
e' 



COS^ : 



b'-(b' k)k 



(A3) 



|e'||b'-(b'-k )k | 

e' b^ 

e"ea [b^-ib^ika/Wka)^] 

e"ba 

e"e„[b^-{b''ka/u»ka)^] 



where we repeatedly used the fact that e" and b" are purely 
space-like in the jet frame. Via explicit computation, we have 
in the lab frame 

e^gfl = e • e = 1 - (z • k) 



e"ba = e-h-^^k^bb 



Therefore, 



b"ka 

u^kb 



cos^ 



e-b+r/3-e^cosi? 
-bcosd . 



e-(b+r/3/7cos?9) 



(A4) 



(A5) 



1 - (z • k) ^sini? 
which reduces to Equation (12) upon noting that e = (l -kk) • z 
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